Quantum fluid effects and parametric instabilities in microcavities 
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We present a description of the non-equilibrium properties of a microcavity polariton fluid, in- 
jected by a nearly-resonant continuous wave pump laser. In the first part, we point out the interplay 
between the peculiar dispersion of the Bogolubov-like polariton excitations and the onset of polari- 
ton parametric instabilities. We show how collective excitation spectra having no counterpart in 
equilibrium systems can be observed by tuning the excitation angle and frequency. In the second 
part, we explain the impact of these collective excitations on the in-plane propagation of the polari- 
ton fluid. We show that the resonant Rayleigh scattering induced by artificial or natural defects is a 
very sensitive tool to show fascinating effects such as polariton superfluidity or polariton Cherenkov 
effect. We present a comprehensive set of predicted far- field and near- field images for the resonant 
Rayleigh scattering emission. 



I. PREFACE 

In the last decade, the research group of Professor 
Marc Ilegems at EPFL has been working intensively and 
enthusiastically on the physics and device applications 
of artificial photonic systems, such as semiconductor mi- 
crocavities and photonic crystals. In this Festschrift pa- 
per, we are going to present a theory of some exotic 
physical properties of coherently excited semiconductor 
microstructures in the strong exciton-photon coupling 
regime. We hope that the rich phenomenology here de- 
scribed will contribute to a very pleasant celebration of 
his 65th birthday. 



II. INTRODUCTION 

The behavior of a quantum fluid has played an impor- 
tant role in many fields of condensed matter and atomic 
physics, ranging from superconductors to Helium flu- 
ids pj and, during the last decade, Bose-Einstein con- 
densates of cold trapped atoms Q- One of the most 
dramatic manifestations of macroscopic coherence of an 
interacting many-body system is superfluidity 

In this paper, we will provide a comprehensive theo- 
retical analysis of the predicted non-equilibrium propa- 
gation properties of a two-dimensional gas of polaritons 
in a semiconductor microcavity in the strong exciton- 
photon coupling regime 0, Thanks to their pho- 
tonic component, polaritons can be coherently excited 
by an applied laser field and detected through the emit- 
ted light. Thanks to their excitonic component, polari- 
tons have strong binary interactions, which have been 
shown to produce spectacular and rich polariton amplifi- 
cation effects through matter-wave stimulated collisions 
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Here, we are going to present a detailed discussion of 
the interplay between the peculiar polariton collective ex- 
citations and the rich variety of parametric instabilities, 
which occur in presence of a nearly-resonant continuous 
wave pump laser. In addition, we are going to discuss 
the impact on the propagation properties of a moving po- 
lariton fluid and analyze in detail the superfluid regime, 
which we predicted in a very recent Letter p7|. As our 
system is a strongly non-equilibrium one, the polariton 
field oscillation frequency is not fixed by any equation 
of state relating the chemical potential to the polariton 
density, but it can be tuned by the frequency of the ex- 
citing pump laser. This opens the possibility of having a 
collective excitation spectrum which has no counterpart 
in usual systems close to thermal equilibrium. In par- 
ticular, we will analyze the propagation in presence of a 
static potential (either for the photonic or exciton com- 
ponent), which is known to produce resonant Rayleigh 
scattering (RRS) of the exciting laser field |lil^l2&l3C|. 
Superfluidity of the polariton fluid manifests itself as a 
dramatic collapse of the resonant Rayleigh scattering in- 
tensity when the flow velocity imprinted by the exciting 
laser beam is slower than the interaction-induced sound 
velocity in the polariton fluid. Furthermore, a dramatic 
reshaping of the RRS pattern due to polariton-polariton 
interactions can be observed in both momentum and real 
space even at higher flow velocities, e.g. with the appear- 
ance of Cherenkov- like patterns. We will present a rich 
set of predicted far-field and near-field images of the res- 
onant Rayleigh scattering emission. 
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III. HAMILTONIAN AND POLARITON 
MEAN-FIELD EQUATIONS 

In order to describe a planar microcavity containing a 
quantum well with an excitonic resonance strongly cou- 
pled to a cavity mode, we will consider the following 
model Hamiltonian 31]: 



n = ( 

ij={X,C} 

+ ^ |dx*] c (x)*t c (x)* x (x)*x(x)+ 

+ J fix hF p (x, t) #^(x) +h.c. , (1) 

where x is the in-plane spatial position and the field op- 
erators \I/jcc(x) respectively describe excitons (X) and 
cavity photons (C). They satisfy Bose commutation 
rules, [#j(x), *t(x')] = c5 2 (x-x') c%. Note that, for sim- 
plicity, we will limit our treatment to the case of polariton 
modes with the same circular polarization, which can be 
excited by a circularly polarized pump. The approach 
here presented can be generalized to the spin-dependent 
case by considering appropriate spin-dependent exciton- 
exciton collisional potentials. 

The single-particle Hamiltonian h° reads 



wc(-iV) 



(2) 



where u>c (k) = \A + k 2 //c 2 is the cavity mode energy 
dispersion as a function of the in-plane wavevector k and 
k z is the quantized photon wavevector along the growth 
direction. is the Rabi frequency of the exciton-cavity 
photon coupling. In the following, we will consider a flat 
exciton dispersion Wx(k) = lux- In this framework of 
coupled harmonic oscillators, polaritons simply arise as 
the eigenstates of the linear Hamiltonian w LP(£/p)(k) 
denotes the dispersion of the lower (upper) polariton 
branch [Fig^a)]. 

The term proportional to F p (x,t) in Eq. repre- 
sents an applied coherent laser pump spot, which drives 
the cavity and injects polaritons. In the following, we 
will consider the case of a monochromatic laser field 
of frequency lo p and plane-wave profile with in-plane 
wave- vector k p . The in-plane wave- vector is linked to 
the incident direction by the simple relationship k p = 
sin#p Wp/c, 9 P being the pump incidence angle with re- 
spect to the growth direction. This means that an oblique 



pump incidence generates a polariton fluid with a non- 
zero flow velocity along the cavity plane. For huj p = 1400 
meV, an in-plane wavevector k p — 1 (fim) -1 corresponds 
to a pump incidence angle of 8.1°. 

The nonlinear interaction term in Eq. is due to 
exciton-exciton collisional interactions and, as usual, is 
modelled by a repulsive (g > 0) contact potential. The 
anharmonic exciton-photon coupling has a negligible ef- 
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FIG. 1: (a): Linear dispersion of the Lower (LP) and Upper 
(UP) Polariton branches as a function of k x (k y = 0). Cavity 
parameters: Tiwx = Tiujc = 1.4 eV, 2TiQ.r — 5meV. (b) 
Zoom of the LP branch. The dashed line depicts the parabolic 
approximation around the bottom of the dispersion. 



feet in the regime considered in the present study [3l| 
and will be neglected for sake of clarity. 

Finally, Vx,c(x) are the single particle potentials act- 
ing on the exciton and photon fields respectively. These 
potentials break the translational symmetry along the 
cavity plane. The exciton potential Vx(x) can be due to 
natural interface or alloy disorder in the semiconductor 
quantum wells due to unavoidable growth imperfections. 
The photonic potential Vc (x) can be due to fluctuations 
of the cavity length or imperfections in the Bragg reflec- 
tors (photonic disorder [241]). More interestingly, Vc(x) 
can be designed and created deliberately by means of 
lithographic techniques. 

Within the mean-field approximation, the time- 
evolution of the mean fields ?/>x,c( x ) = (^JC,c( x )) under 
the Hamiltonian is given by: 



d ( V>x(x) 



dt V V>c( x ) 



Ulpt) 



Vx(x) 



7|<Mx)| 5 






Vc (x) - i7c 



i>c ( x ) 



(3) 



Using the language of the quantum fluid commu- nity, these are the Gross-Pitaevskii equations 0] for our 
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FIG. 2: Exciton mean-field energy Tig^x | 2 (meV) as a func- 
tion of the incident pump intensity (arb. units). Cavity pa- 
rameters: h-yc = hjx = 0.1 meV, Tilox = — 1.4eV, 
2 TiQ.r — 5meV. (a) Bistability curve obtained with the exci- 
tation parameters: k p = 0.314 /zm~ (well in the parabolic re- 
gion near the bottom of the LP dispersion), hu! p — hiJLp(k P ) ~ 
0.47 meV. Circles depict the calculated stable points, while 
the dashed line represents the unstable branch. The thresh- 
old points A and B are, respectively, due to a single-mode 
(Kerr) or a multi-mode (parametric) instability, (b) Op- 
tical limiter curve obtained with the same k p , but with 
hcjp — hu)Lp(k. P ) ~ — 0.47meV. In this case, all stationary 
solutions are stable. 



are the non-equilibrium analogous of the state equation, 
which in equilibrium systems links the chemical potential 
to the particle density. Importantly, we stress that while 
the oscillation frequency of the condensate wavefunction 
in an isolated gas is equal to the chemical potential fi 
and therefore it is fixed by the equation of state, in the 
present driven-dissipative system it is equal to the fre- 
quency uj p of the driving laser and therefore it is an ex- 
perimentally tunable parameter. Hence, the microcavity 
polariton system allows us to explore a regime, which is 
not accessible in systems close to thermal equilibrium, 
such as the ultracold trapped atoms. 



V. LINEARIZED BOGOLIUBOV-LIKE THEORY 

As usual in the theory of nonlinear systems, stability 
of the solutions of Eqs. 1)415(1 with respect to fluctua- 
tions has to be checked by linearizing Eq. |J3J around 
the stationary state. Perturbations can be produced by 
classical fluctuations of the pump field, quantum noise 
of the exciton and photon fields as well as the presence 
of perturbing potentials Vc,x(x), which have not been 
considered by the plane- wave solutions in Eqs. (11151 . 

In the stability region, the linearized response of the 
system to a weak perturbation is analogous to the cele- 
brated Bogoliubov theory of the weakly interacting Bose 
gas 0. Let us define the slowly varying fields with re- 
spect to the pump frequency as 



cavity-polariton system. The quantities jx and 7c rep- 
resent the homogoneous broadening of the exciton and 
photon modes respectively. In the present work, we will 
be concerned with an excitation close to the bottom of 
the LP dispersion, i.e. the region most protected |5j from 
the exciton reservoir, which may be otherwise responsible 
for excitation- induced decoherence |32| ■ 



IV. STATIONARY SOLUTIONS IN THE 
HOMOGENEOUS CASE 

In the homogeneous case (i.e., Vx,c( x ) = and trans- 
lational invariance along the plane), we can look for 
spatially homogeneous stationary states of the system 
in which the field has the same plane-wave structure 
i/jx.c( x , t) = exp[i(k p x — ui p t)] ipx c as the incident laser 
pump field. The mean-field equations 

(cjx(k) -u P -^x+g \rx\ 2 ) n + n R rj = (4) 

(wc(k P ) -u p - l -ic) 1>c + Mr 4>x = -F P , (5) 



<5</>i(x, t) = Sipi(x, t) exp(iu p t) , (6) 

and let us consider the four-component displacement vec- 
tor 

^(x,i) = (^x(x,i),^c(x,i),^(x,0,^c(x,i)) T • 

The equation of motion for the four-component displace- 
ment vector reads 

= £ ■ 5$ + fpert, (8) 

where f per t is the inhomogeneous source term produced 
by the perturbation. The expression for f per t depends on 
which kind of perturbation is considered and will be given 
explicitly later for the case of a perturbation induced by 
the single particle potentials Vc7,x( x )- The linear opera- 
tor C is 
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FIG. 3: Same plot as in Fig. |5J but with pump in-plane 
wavevector k v = 1.5 /im -1 (close to the inflection point of 
the LP dispersion). Pump frequency: (a) hui p — huiLp(k p ) = 
0.47 meV, (b) hui p — huLp(k P ) = — 0.47 meV. In contrast 
to Fig. |5] here both threshold points C and D are due to 
parametric instabilities. 



polariton energy (u> p > uir.pCkn))- In this c ase, there is a 
clear S-shaped bistability curve El ES E5- The un- 
stable branch, determined through the eigenvalues of the 
linear operator C, have been depicted with a dashed line, 
while the stable points are represented by circles. The 
threshold points A and B are due to a Kerr and to para- 
metric instability respectively. By comparison, in Fig. 
|2Ib), we have shown the same quantity, but for a red- 
detuned laser frequency (ui p < ULp(k p )). In this case, 
the polariton system behaves as an optical limiter[3^|. 
the absorption is highly sublinear and all the points are 
stable. To give a more complete picture, we report in Fig. 
[31 the analogous calculations, but with a larger wavevec- 
tor (k p = 1.5 /im -1 ), close to the inflection point of the 
LP dispersion. It is apparent that, while the shape is 
analogous, the boundary between the stable and unsta- 
ble branches is modified. In particular, the threshold 
points C and D are both due to parametric instabilities. 
Note that nice hysteresis loops due to polariton bistabil- 
ity have been recently experimentally demonstrated in 
the case k p = [34| and for a pump wavevector close to 
the inflection point of the LP dispersion 36] . 

B. Complex energy of the collective excitations 



A. Stability of the stationary solutions 

The stability of the solutions of Eqs. (1 1151 can be de- 
termined by calculating the imaginary parts of the eigen- 
values of the operator L. If all the eigenvalues have 
negative imaginary parts (i.e., as it happens in the non- 
interacting case), then the solutions are stable. Other- 
wise, an instability occurs. If the polariton instability in- 
volves only the pump mode, we have the analogous of a 
Kerr instability. If the instability is due to pairs of modes 
(formation of the so-called signal-idler pairs), we have a 
parametric instability (in the field of quantum fluids, this 
kind of dynamical instabilities are generally known as 
modulational instabilities |33| ) . In Fig. we have plot- 
ted the stationary solutions for the exciton mean-field en- 
ergy hg \ipx 1 2 (meV) as a function of the incident pump 
intensity (arb. units) for realistic microcavity parameters 
and with a small pump wave- vector (k p = 0.314 /im -1 ), 
close to the bottom of the LP dispersion. In Fig. Ufa), 
we have considered the case of a pump frequency, which 
is blue-detuned with respect to the unperturbed lower 



The spectrum of the collective excitations (Bogoliubov 
modes in the quantum fluid terminology) can be obtained 
from the eigenvalues of the operator C As the system is 
translationally invariant along the plane (we are consider- 
ing the homogeneous case Vx = Vc = 0) , the wavevector 
k is a good quantum number and the eigenvectors of C 
have a plane-wave form 
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ikx 



i(k— 2k p )x 
= i(k— 2k p )x 



satisfying the reduced eigenvalue equation 



(wf(k)- Wp )l-£(k ) k p ) 



where 



(10) 



0, (11) 
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For each k, the spectrum is composed by four branches. 
For each polariton branch j £ {LP, UP}, two ± branches 
exist, which are related by the symmetry 

wjQt) = 2u p - w/(2kp - k). (13) 

Now, we wish to point out and list clearly the relevant 
properties and symmetries in this problem. These prop- 
erties will be later discussed in detail through a set of 
comprehensive examples and elucidations. 

(i) The collective excitations are characterized by the 
pump-induced coherent coupling between a generic mode 
with wavevector k and the "idler" wavevector 2k p — k. 
This corresponds to the elementary process {k p ,k p } — > 
{k, 2k p — k}, i.e. the conversion of two pump excitations 
into a signal-idler pair (to use the quantum optics ter- 
minology of parametric oscillators) with the same total 
momentum. 

(ii) The "idler" branch w£p(k) is the "image" of 
the ordinary branch u>^ P (k) under the simultaneous 
transformations 38] k — > 2k p — k and uj — > 2uj p — u>. The 



same relationship holds for the UP branches wjjp(k). 

(hi) The matrix £(k, k p ) is characterized by an anti- 
hermitian coupling between k and 2k p — k. This feature is 
typical of parametric wave-mixing coupling (in the quan- 
tum optics language) or Bogoliubov theory (using the 
quantum fluid literature terminology). 

(iv) The four branches of eigenvalues u>^j P LP (k) are 
complex. As the real parts, the imaginary parts of the 
eigenvalues depend both on k and on the pump param- 
eters. In the stability region, all imaginary parts are 
negative. 

(v) In case of resonant excitation of the lower branch, 
provided that the interaction energy gIV'xl 2 ^ s mucn 
smaller than the polaritonic splitting u>up~u>lp, there is 
no significant mixing between the LP and UP branches. 
Hence, a simplified approach consists in neglecting the 
contribution of the upper branch. With this approxi- 
mation, the branches w^p(k) are the eigenvalues of the 
simplified matrix 



£ip(k,k p 



w L p(k) + 2g LP \<ip 



LP\ 

-9lpVl s p 2 



t7Lp(k) 
2 



9 LP Vlp 

2uj p - {u LP (2k p - k) + 2g LP |Vfp| 



»7Lp(k) 
2 



(14) 



where the stationary lower polariton field is written as 
a linear superposition of the exciton and cavity photon 
fields, namely 

rffp = X LP (k p )r x s + C L p(k p )r c s , (15) 

being | Alp (k p )| 2 and |C,Lp(kp)| 2 the exciton and pho- 
ton fractions of the lower polariton mode with the pump 
wavevector. The effective interaction strength 

g LP = g \X L p(k p )\ 2 X LP (k)X LP (2k p - k) (16) 

takes into account for the exciton fraction of the involved 
lower polariton modes (pump, signal and idler). The 
polariton linewidth is given by 7Lp(k) = \XLp(k)\ 2 "fx + 
|C L p(k)| 2 7c . 

(vi) When the diagonal elements of C^pik, k p ) are 
equal, it is easy to verify that 9?[a/£p(k)] = ^t[uj^ p (k)], 
while Sy[o;Jp(k)] =^ Q[w£ p (k)]. This means that the para- 
metric coupling produces a splitting of the imaginary 



parts of the two LP branches, while the real parts are the 
same. If the difference between the diagonal elements of 
£Lp(k, k p ) is small compared to the coupling g^p IV'lpI 2 * 
then the same property holds. In other words, the dis- 
persions of the two branches uj~^ P (k) and ui^ p (k) stick to- 
gether, while their imaginary parts are split. One branch 
is narrowed with respect to the linear regime, while the 
other is overdamped. Analogous properties occur for the 
exact eigenvalues of the 4x4 matrix £(k, k p ), which are 
reported in all the figures of this paper. 



1. Excitation near the inflection point of the LP dispersion 

In the following, we will show the exact eigenval- 
ues (obtained by numerical calculations) of the matrix 
£(k, k p ) in Eq. I|12fl a s a function of the excitation pa- 
rameters (pump frequency, wavevector and intensity). 
We will focus on the subtle interplay between the dra- 
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FIG. 4: (a) Exact energy dispersions ^t[hu}^ p up ] of the four 
polariton Bogoliubov branches measured with respect to the 
pump photon energy TvjJ v (meV). (b) Corresponding imagi- 
nary parts (meV). Note that negative imaginary parts imply 
stability. Excitation parameters: k p — 1.5 /im" 1 (along the 
x-axis), TitOp - huj L pCkp) = 0.107 meV, hg^xl 
Cavity parameters as in the previous figures. 



M*xl 2 = 0.05meV. 
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(a) Exact energy dispersions ^St[hu}^ 
polariton Bogoliubov branches measured with respect to the 
pump photon energy hu> p (meV). (b) Corresponding imagi- 
nary parts (meV). Excitation parameters: k p — 0.314 fjxn~ 
(along the x-axis), fkj p — TiuiLp(k- p ) — 0.47 meV, hg^x | 2 = 
1.02meV. Note that this case is the precursor of a Kerr in- 
stability, because the imaginary parts are modified at the 
pumped mode only. This stable point is close to the inversion 
point A in Fig. |5Ja) . 
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FIG. 5: (a) Exact energy dispersions ^[hujj^p] of the lower 
polariton Bogoliubov branches measured with respect to the 
pump photon energy hui p (meV). (b) Corresponding imagi- 
nary parts (meV). Excitation parameters: k p = 1.5 jitm - 
(along the x-axis), huj p — hu)Lp(it p ) = 0.47 meV, fig^x] 2 = 
0.699meV. Note that here the stationary solution is unstable, 
because there are modes with positive imaginary parts. This 
unstable point is close to the point C in Fig. EI a). 



matic modification of the energy dispersions (depend- 
ing on the real part of the eigenvalues) and the onset of 



the parametric instabilities (depending on the imaginary 
part). 

As a first example, we consider the case of nearly- 
resonant excitation close to the inflection point of the 
LP dispersion (see Fig. The pump frequency has 

been taken slightly blue-detuned with respect to the po- 
lariton energy in the linear regime. In Fig. 0Ja) the exact 
dispersions of the four polariton Bogoliubov branches is 
shown. The upper polariton branches are energetically 
far away and play a negligible role, while the relevant 
physics concerns the lower polariton branches ui^p only. 
The corresponding imaginary parts are shown in lljTb) • 
It is apparent that there is a dramatic modification of 
the imaginary part around the wave-vectors k x = and 
k x = 2k p . Although the stationary solutions are here 
stable (negative imaginary parts), we can see that we are 
close to a parametric instability. In fact, there is one 
branch, whose imaginary part is not far from zero. Note 
that the imaginary parts are split at k p as well, even if 
in a much weaker way. This is a precursor of a Kerr (or 
single-mode) instability. 

In Fig. [S] we give another example, with the same 
excitation parameters as in Fig. 131(a) and with an exciton 
mean-field energy hg\^ s ^\ 2 — 0.699meV. In this case, we 
have an unstable solution, because, as shown in FiglSfb), 
there are modes with positive imaginary parts. On the 
bistable curve of Fig. EJa), we are here just beyond C. 
The instability being of parametric (or multi-mode) type, 
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FIG. 7: (a) Exact energy dispersions ^[huj^p] of the lower 
polariton Bogoliubov branches measured with respect to the 
pump photon energy hui p (meV). (b) Corresponding imag- 
inary parts (meV). Excitation parameters: k p — 0, hui p — 
fuuLPikp) = 0.532 meV, hg\$x\ 2 = 1.2meV. 
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FIG. 8: (a) Exact energy dispersions ^[hujj^p] of the lower 
polariton Bogoliubov branches measured with respect to the 
pump photon energy hui p (meV). (b) Corresponding imagi- 
nary parts (meV). Excitation parameters: k p — 0.314 /zm~ 
(along the x-axis), Tllj p ~ hcjLp(k. P ) = 0.04meV, ftc/l^xl 2 = 
0.04meV. 



note that the point C is close but does not coincide with 
the inversion point of the bistable curve. In Fig. |Sfa), 



we can see that the branches 



J LP 



stick together in the 



wavevector region where the parametric instability takes 
place. 



2. Excitation near the bottom of the LP dispersion 

Here, we consider the case of a smaller pump excitation 
wavevector and energy , such as to excite the LP branch 
close to the bottom of its dispersion. In this region, as 
shown by Fig. 0b), the dispersion of the unperturbed 
lower polariton branch is parabolic. In order to stress the 
non-trivial effects here predicted, we start by showing a 
spectacular case, depicted in Fig|5] The excitation pa- 
rameters correspond to Fig. |2Ja) with the exciton mean- 
field energy hg\^x\ 2 — l-02meV, i.e., a point close to 
the threshold point A for the Kerr instability. The dis- 
persions of the branches uj^p in FigEt a ) have a corner 
at the pump wavevector. This dispersion is reminiscent 
of the celebrated Bogoliubov linear dispersion in super- 
fluid helium and in the atomic condensates. If we look 
at FigHJb), we realize that the modification of the imag- 
inary part is peaked around the pump wavevector itself. 
Using the quantum optics language, this is the precur- 
sor of a Kerr instability, because it involves the pumped 
mode only. 

In Fig. [7] we give another example, which has no 
analog in equilibrium systems. Here, we have a blue- 
detuned pump at normal incidence, namely k p = 0, 



huj p - hw L p(k p ) = 0.53 meV and Tig\^f \' A = 1.2meV. 



In Fig. 0a), we can clearly see that the dispersion of the 
polariton collective excitations is flat around the pump 
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FIG. 9: (a) Exact energy dispersions ^H[huij^ p } of the lower 
polariton Bogoliubov branches measured with respect to the 
pump photon energy hco p (meV). (b) Corresponding imagi- 
nary parts (meV). Excitation parameters: k p — 0.314 fj,m~ 
(along the x-axis), hiv p — hLu L p(k. p ) = 0.47 meV, hg\^x\ 2 = 
0.04meV. 



wave- vector. For the parameters of Fig. Efb), the imag- 
inary parts are all negative, which implies stability, but, 
as in FigHJb), we are not far from the onset of a Kerr 
instability. Note that in Fig. |SJ we have an analogous 
situation, but with a finite pump wavevector. As shown 
by Fig. |HIa), the branches stick together around the 
pump wavevector, with a dispersion exactly linear. 
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In Fig. HJJa), we show the dispersions for a stable point, 
which is close to the threshold point B in Fig. Ufa). 
Note that here the exciton mean-field energy is consider- 
ably smaller than the pump detuning. Hence, the branch 
sticking occurs in a limited portion of momentum space, 
where the imaginary parts are affected, as shown in Fig. 
MP). 

Finally, in Fig. HUf a) . we give another different ex- 
ample, with a full gap between the branch lo\ p and the 
branch oJ^p- In Fig. llOf b) . we can see that the imagi- 
nary parts of the eigenvalues are unchanged with respect 
to the linear regime. Indeed, for these excitation param- 
eters (see the caption of Fig. 1100 . the microcavity has an 
optical limiter behavior similar to the one in Fig. HJb) 
and no instability occurs if the pump intensity is further 
increased. 



variety of interaction-induced polariton instabilities, we 
present here a simple approximated approach, which al- 
lows us to grasp effectively the physics contained by the 
eigenvalues of the matrix in Eq. I|12f) . I n particular, 
we consider the case of negligible mixing with the UP 
branches, which allows us to focus our analysis on the 
simpler 2x2 matrix in Eq. (|14|) instead of the 4x4 
matrix in Eq. (|12fl . Moreover, we consider a pump ex- 
citation close to the bottom of the LP dispersion, where 
the dispersion is approximately parabolic (see Fig. db)), 
i.e., 



^Lp(k) ~ lolp{0) 



hk 2 

2m LP 



(17) 



3. Simplified analytical model for excitation close to the 
bottom of the LP dispersion 

Now, after having shown a few examples of the rich 
spectra of non-equilibrium collective excitations and the 



where vtllp is the effective mass of the LP dispersion. 
Under these assumptions, the spectrum of the LP Bo- 
goliubov excitations can be approximated by the simple 
expression 



J 



J LP 



5k- 



IJLP 



± 



sj (2 g L p\ip S Lp\ 2 + VSk - Ap)(77 Sk - A p ), 



(18) 



where Sk — k — k 



■v 



V6k 



hSk 2 

2m L p 



(19) 



the pump mode flow velocity is v p = Tikp/rriLP and the 
interaction-renormalized pump detuning 



■u LP {k p ) - 9lp\Vlp\ 



(20) 



In the case of resonant exciton-photon coupling (i.e. 
w c(0) = u>x) and small wavevectors, the excitonic 
fraction of the lower polariton mode is approximately 
0.5. Under these assumptions, we have I^lpI 2 ~ 
2\ip^\ 2 and g LP » g/A [see Eqs. JTSITffl l. Therefore, 
gLp\ip S Lp\ 2 ~ 0.5 g\ipx\ 2 ' i- e - the mean-field interac- 
tion energy "felt" by the lower polariton is half of the 
mean-field energy for the exciton field. Note that there 
are three different cases, according to the value of A p 
defined in Eq. lt5Ul). 

(i) Ap = 0. In this resonant situation, the ± branches 
touch at k = k p . The effect of the finite flow velocity 
Vp is to tilt the standard Bogoliubov dispersion via 
the term 5k ■ v p . While in the non- interacting case the 
dispersion is parabolic, in the presence of interactions 
[FigUJa)] its slope has a discontinuity at k = k p . On 
each side of the corner, the + branch starts linearly with 



group velocities respectively given by v 
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FIG. 10: (a) Exact energy dispersions ^t[hojj^ p ] of the lower 
polariton Bogoliubov branches measured with respect to the 
pump photon energy hcu p (meV). (b) Corresponding imagi- 
nary parts (meV). Excitation parameters: k p — 0.314 /zm~ 
(along the x-axis), huj p — hoJLpikp) = 0.25 meV, hg\^x | 2 = 
0.6meV. 



being the usual sound velocity of the interacting Bose gas 
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ft 9Lp\^Lp\ 2 /m L p . 



(21) 



On the hysteresis curve of Fig^a), the condition A p = 
corresponds to the inversion point A. 

(ii) A p > 0. In this case, the argument of the square 
root in (|18fl is negative for the wavevectors k such that 
A p > rjsk > Ap— 2 g\^x I 2 - ^ n this region, the ± branches 
stick together |3l| (i.e. 3?[w^ P ] = 5R[w^ P ]) and have an 
exactly linear dispersion of slope v p as in Fig. OJa) and 
in Fig. Qa). The imaginary parts are instead split, with 
one branch being narrowed and the other broadened [3ll 
l38| . Increasing further the pump density, the multi-mode 
parametric instability |38| sets in, corresponding to the 
point B in Fig. Ufa) ■ 

(hi) Ap < 0. In this case, as it is shown in Fig llOl 
the branches no longer touch each other at k p and a full 
gap between them opens up for sufficiently large values 
of |Ap|. In Fig. EJa), the region Ap < is indicated. 



sponse to a perturbation. In particular, we shall con- 
sider here a moderate static disorder as described by the 
potential Vc,x(x). In this case the perturbation source 
term for the equations of the linearized theory (see Eq. 
IHJ) is the time-independent quantity 



Vb(x)#? 
-Vx(x) 



ss * 
X 
s * 



(22) 



The induced perturbation of the exciton and photon 
fields is given by the expression 



fc(x) = -C- 1 ■ / d (x) 



(23) 



VI. RESPONSE TO A STATIC POTENTIAL: 
RESONANT RAYLEIGH SCATTERING 

The dispersion of the polariton elementary excitations 
is the starting point for a study of the microcavity re- 



We remind you that, as shown by Eq. £ is an 

operator depending on the two-dimensional spatial gra- 
dient V. It is convenient to perform a spatial Fourier 
transform, which leads to the algebraic result 



^c(k) 
l* x (2k p -k) 
l* c (2k p -k) / 



A" 



-(£(k, kp) - huj p 



V x (k) 
Vb(k) 4% 

-y x (k-2k p ) 

-V x (k-2k p ) 



ISS * 

D X 
tss * 



(24) 



where the eigenvalues of the matrix £(k, k p ), defined in 
Eq. (|12fl , are the 4 branches of polariton Bogoliubov ex- 
citations. The perturbation potentials Vc,x(x) break the 
planar translational symmetry of the microcavity system, 
thus exciting polariton modes with in-plane wavevectors 
different from the pump wavevector k p . However, being 
Vc,x(x) static, the resonant excitation concerns only Bo- 
goliubov modes whose frequency is equal to uj p (within 
the polariton homogeneous linewidth). The observable 
quantity is 



I RRS (k) oc |^ c (k)| 2 , (25) 



i.e., the perturbation-induced intensity of the photonic 
field, which is proportional to the far-field images of the 
resonant Rayleigh scattering signal 0, 0> HH ■ 



A. Weak excitation regime and elastic RRS ring 

In the following, we will show a few applications of Eq. 
(|2~p . using the perturbation potentials depicted in Fig. 
mi Here, we have considered a single photonic defect 
(depth 1 meV, width 1.5/im) for the in-planc photonic 
potential Vc(x) [see Fig lllf a'l] and a disordered excitonic 
potential [see Fig lllf b - )]. We point out that these poten- 
tials are just an example and that Eq. I|24|l can be read- 
ily applied for an arbitrary set of perturbation potentials 
(see, e.g., the cover of this special volume). In the nu- 
merical applications here reported, we have considered 
a 256 x 256 spatial grid and a squared box (400/im x 
400/im), with the photonic dot at the center of the box. 

In Fig. El we show the results for the weak excita- 
tion regime, where the many-body effects produced by 
polariton-polariton interactions are negligible. In this 
linear regime, we have the conventional unperturbed dis- 
persions of the lower and upper polariton branches. In 
Fig- EI a ) j we have considered the case of resonant exci- 
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FIG. 11: Single particle potentials considered in the numerical 
calculations. Note that this is a zoom around the origin of 
the 400>m x 400/Um box (with a 256 x 256 grid) used for the 
numerical calculations. Top panel: Photonic potential Vc(x) 
(meV) (it can model an artificial or natural point defect at 
the origin x = y = 0). Bottom panel: Excitonic potential 
Vx (x) (disordered spatial fluctuations of the exciton energy) . 
The gray color scale is different with respect to the top panel. 
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FIG. 12: Energy dispersion of the LP branches in the 
weak excitation regime. (b) Intensity (arb. units) of the pho- 
tonic resonant Rayleigh scattering signal |<50c(k)| 2 . Ex- 



citation parameters: k p 



0.4 fim 

ss|2 



(along the x-axis) 



Tiu p - huj L p(k p ) = OmeV, = O.OOOlmeV. The pho- 

tonic and exciton potentials are those shown in Fig 1111 The 
speckles of the elastic RRS ring are due to the disordered ex- 
citonic potential. The chosen plot range is such that a few 
speckles saturate the gray scale. 



Normalized photon intensity 




tation close to the bottom of the LP branch, where the 
dispersion is parabolic. In Fig. I12f b). the intensity of 
the resonant Rayleigh scattering is shown, displaying the 
well known elastic ring. In fact, in the linear regime, the 
solutions of the equations u>Lp(k) = lu p are k-points on 
a circle, because the unperturbed polariton dispersion 
depends on |k| only. The speckles on top of the elas- 
tic ring are due to the random nature of the excitonic 
potential. The width of the ring is due to the finite homo- 
geneous broadening of the polariton modes. Note that, 
in order to excite the elastic ring corresponding to the 
wavevector k p , the Fourier component Vc,x0^p) of the 
static potentials need to be finite. This condition is eas- 
ily fulfilled by a typical excitonic disordered potential or 
by a photonic defect whose width is of the order of 1 /im. 
In Fig. 1131 we show the corresponding spatial pattern. 
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FIG. 13: Spatial profile of the normalized cavity photon den- 
sity, i.e., Ic{~x) / Ic° m ■ Excitation parameters and potentials 
as in Fig. 1121 The coherent diffusion pattern induced by the 
point defect at the origin is the main feature on top of the 
random landscape produced by the exciton disorder. 
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FIG. 14: Superfluid regime. Same parameters as in Fig. 1121 
but with fiu p -?uj L p(kp) = 0.467 meV, %|*x| 2 = lmeV. In 
this superfluid regime, the RRS elastic ring has collapsed. 
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FIG. 15: Spatial profile of the normalized cavity photon den- 
sity in the superfluid regime. Excitation parameters as in Fig. 
1141 To compare with the normal (weak) excitation regime, see 
Fig. H31 

Precisely, we have plotted the normalized quantity 
^c(x) |# 7 V k »-* + ^c(x)| 2 

i.e., the total photon field intensity (homogeneous solu- 
tion + potential-induced perturbation) normalized to the 
intensity of the homogeneous solution without the poten- 



tial. For the considered potentials in Fig. 1111 the dom- 
inant feature is due to the photonic point defect. The 
polariton plane wave driven by the pump is coherently 
scattered by the photonic defect (located at the posi- 
tion x = y = 0), producing a peculiar interference pat- 
tern, characterized by parabolic wavefronts. In fact, the 
polariton field scattered by the point defect is a cylin- 
drical wave. Hence, if we consider only one defect and 
k p is along the x-direction, the total field has the form 

/(x) = e lkpX + {3e zkp V /x2+y2 . The constant phase curves 
are given by the condition k p x + k p yjx 2 + y 2 — 2im, 
whose solutions describe parabolic wavefronts with a 
symmetry axis oriented along the x-direction, as nicely 
depicted by the exact solution in Fig. Due to the 
presence of the exciton potential, additional disordered 
features are superimposed on the main interference pat- 
tern produced by the photonic point defect. 



B. Superfluid regime 

In presence of interactions, we have seen that the spec- 
trum of polariton Bogoliubov excitations is dramatically 
different from the unperturbed case. This manifestation 
of polariton many-body physics can be probed in a sen- 
sitive way by the resonant Rayleigh scattering emission. 
In Fig. 1141 we start by considering the most spectacu- 
lar regime of polariton superfluidity. This regime can be 
achieved when the pump is resonant with the interaction- 
renormalized polariton dispersion at the pump wavevec- 
tor (A p = 0) and when the sound velocity c s [see Eq. 
I)2ip] of the interacting polariton fluid is larger than the 
flow velocity v p = Tik p /rriLP imprinted by the pump 
beam. This situation is more favorable to obtain for ex- 
citation close to the bottom of the LP dispersion, imply- 
ing smaller pump flow velocity v p and smaller excitation 
density necessary to have c s > v p . As depicted by Fig. 
H"4T a1. in this case, the equation 5R[o;^p(k)] — lu p = has 
no solutions for k ^ k p , meaning that no final states are 
available for the elastic scattering induced by the static 
potential. As a dramatic consequence, the elastic ring 
in Fig. I12r b) collapses. As shown in Fig. I14f b). only 
a weak emission around the pump wavevector k p is left, 
due to non-resonant processes, which are allowed by the 
finite broadening of the polariton modes. 

The real space pattern is shown in Fig. 1151 showing 
that the effect of the disorder remains localized around 
the defect positions. Hence, the polaritonic propagation 
is superfluid. In analogous way to liquid Helium and 
atomic condensates [HQ, we can state that the polariton 
fluid has a superfluid behaviour according to the Landau 
criterion, if and only if both following conditions are sat- 
isfied 

(a) oj^ PUP (k) > uj p for every k ^ k p . 

(b) Lo p > o->lp(0), i.e. there is an elastic ring in the 
weak excitation regime. 

We point out that the condition (b) is necessary to 



12 



have a meaningful definition of polariton superfluidity. 
In fact, if Up < Wj,p(0), already in the weak excitation 
regime there are no real states at the pump energy and 
there is no resonant Rayleigh scattering clastic ring. Note 
that the conditions (a) and (b) are achieved not only in 
the resonant case A p = 0. Within the parabolic approx- 
imation in Eq. I|18|) . conditions (a) and (b) are satisfied 



when A p < and 



-gLp\ip 



LP I 



m LP v^ 
2H 



< A p < g L p\i> 



m LP v v 



LP\ 



(27) 

The calculations here reported show the robustness of 
the superfluid flow with respect to elastic processes, such 
as the scattering on static defects. As the Landau cri- 
terion for superfluidity involves also inelastic processes 
(e.g., emission of crystal phonons or heating of residual 
free carriers), it is important to note that whenever the 
superfluidity condition (a) for elastic scattering is ful- 
filled, then it is satisfied a fortiori also for the inelastic 
channels. In fact, the stability of the mean-field solu- 
tion implies that almost no Bogoliubov quasi-particles 
are present above a stable mean-field solution, so when- 
ever oj^p(k) > uj p , no final states are available for the po- 
lariton system to lower its energy by transferring energy 
to its environment. This is here assumed to be at almost 
zero temperature, so that it can only absorb energy from 
the polariton system. Friction is therefore absent in this 
regime. 

On the other hand, at a finite temperature of the crys- 
tal a thermally excited "normal" component can appear 
in the polariton fluid because of its heating due the inter- 
action with the phonon bath and the residual free carri- 
ers. As it happens in liquid Helium, this normal compo- 
nent experiences a finite friction when flowing onto the 
defect. However, its magnitude vanishes at low temper- 
atures, and in any case it does not affect the superfluid- 
ity of the co-existing superfluid component in a two-fluid 
picture. 



C. Precursors of parametric instabilities and 
branch sticking 

In the case A p > (i.e., pump frequency higher than 
the renormalized frequency of the pumped mode), the 
resonant Rayleigh scattering response is completely dif- 
ferent, as shown in Fig. 1161 and ITS1 In this regime, the 
two LP branches stick together, while there is a splitting 
of their imaginary parts [see, e.g., the analogous situa- 
tion in Fig. ISJand Fig. |SJ. Such a scenario represents the 
precursor of a multi-mode parametric instability, which 
can be triggered by further increasing the excitation den- 
sity. In contrast to the superfluid regime, a deformed 
RRS ring is apparent in Fig. E^b) and Fig. Elb). In 
particular, depending on the topology of the k-space re- 
gion where the LP branches stick, the RRS intensity is 
strongly amplified either on a segment parallel to y in- 
cluding the point k p [Fig lltjt b)] , or around two points of 
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FIG. 16: Same parameters as in Fig. 1121 but with hu> p — 
hcj LP (kp) = 0.1 meV, %|*x| 2 = 0.07meV. Branch sticking 
and amplified RRS are precursors of a parametric instability. 



Normalized photon intensity 




FIG. 17: Spatial profile of the normalized cavity photon den- 
sity. Excitation parameters as in Fig. 1161 Note that the gray 
scale of this plot is saturated in the region around the point 
defect at the origin. 



the straight line parallel to y and passing through the 
point k p [Fig. 118( b)] . The main consequence of this in 
the real-space pattern of Fig. El is an overall amplifica- 
tion of the density modulation induced by the defect, in 
stark contrast with the superfluid regime. In particular, 
note the long " shadow" in the downstream direction with 
respect to the central defect, which extends to relatively 
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FIG. 18: Excitation parameters: k p — 0.314 /im _1 (along the 
x-axis), huj p - huj LP (]s. p ) = 0.47 meV, Tig\^x\ 2 = 0.075meV. 
Note that the gray scale of this plot is saturated around the 
two brightest points. Here, we consider the situation of a 
dominant photonic potential (Vx = 0). 
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FIG. 20: Cherenkov regime. Parameters: u>x = ^c(0) — 
lmeV, k p — 0.7/im -1 (along the x-axis), Tiui p — 7kd£,p(k p ) = 
0.599 meV, %|*f] 2 = lmeV. Here, we consider the situation 
of a dominant photonic potential (Vx = 0). 
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FIG. 19: Spatial profile of the normalized cavity photon den- 
sity. Same parameters as in Fig. ^| 
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FIG. 21: Spatial profile of the normalized cavity photon den- 
sity. Same parameters as in Fig. 1201 The photonic point 
defect produces Cherenkov-like wavefronts. Note that, in or- 
der to show the peculiar shape of the wavefronts, the gray 
scale of this plot is saturated in the region around the point 
defect . 



far distances thanks to the linewidth narrowing effect. 
In Fig^l the shadow of the defect is even more pecu- 
liar, showing a series of fringes parallel to the x direction. 
These can be explained in terms of the interference be- 
tween the pump and the two peaks in k-space shown in 

Figdb). 



D. Cherenkov regime 

Here, we consider the opposite case A p < 0, but with 
an excitation density which is not high enough to enter 
the superfluid regime characterized by the condition in 
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FIG. 22: Parameters: ojx = wc(0), fc p = 0.7/im 1 (along the 
x-axis), hojp - ftwLp(kp) = 0.3 meV, %|*3f | 2 = 0.6meV. 



FIG. 24: Parameters: ujx = u>c(0), k p — 0.7/im 1 (along the 
x-axis), hojp - huj L p(k p ) = 0.2 meV, Tig\^x\ 2 = 0.6meV. 
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FIG. 23: Spatial profile of the normalized cavity photon den- 
sity. Same parameters as in Fig. 1221 The gray scale of this 
plot is saturated in the region around the point defect. 



Eq. i|27|) . For the sake of clarity, as we have already done 
for Fig. 1181 and Fig. 1191 we take here Vx — 0, so to con- 
centrate on the effect of a single defect. This situation 
can be realistically realized, e.g., when there is a single 
photonic defect (natural or artificial) which is dominant 
over the background excitonic disorder. In Fig. I20f a). 
we consider the resonant case (A p = 0), with a polari- 
ton sound speed c s < v p . As shown in Fig. !2UT b) the 
weak excitation elastic RRS ring is replaced by an asym- 



Normalized photon intensity 




-50 -40 -30 -20 -10 10 31 30 40 50 



x [ ji m ] 

FIG. 25: Spatial profile of the normalized cavity photon den- 
sity. Same parameters as in Fig. 1241 The gray scale of this 
plot is saturated in the region around the point defect. 



metric pattern, which is strongly deformed and shows a 
singularity at the pump wavevector. The aperture angle 
26 of the singularity at k p satisfies the simple condition 
cos 9 = c s /vp. In this c s < v p regime where the polariton 
fluid is moving at a supersonic speed, the defect produces 
a peculiar real-space pattern (Fig . I21fl showing linear 
Cherenkov-like wavefronts [39l |40( . The aperture 2<p of 
the Cherenkov angle has the usual value sin<^ = c s /v p . 
This behavior is easily understood from a physical stand- 
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point as follows: a moving fluid propagating along the 
positive x direction in the presence of a static defect is 
equivalent, under a Galilean transformation, to a defect 
moving in the negative x direction in a fluid at rest. This 
situation is a familiar one, as it corresponds to the wave- 
fronts created by a moving duck on the surface of a lake. 
The rounded region on the left-hand side of the k-space 
pattern in Fig. I20r b) was not present in the standard the- 
ory of Cherenkov emission in non-dispersive media |39| 
and it is due to the fact that the Bogoliubov dispersion is 
linear only in the neighborhood of k. p and then it bends 
upwards. A remarkable consequence of this property is 
the oscillatory perturbation shown by the real-space pat- 
tern upstream with respect to the defect, as apparent in 
FigET] 

In the case A p < 0, the branches ui^p an( i ^lp are 
disconnected, as well as the corresponding RRS emission 
pattern, as depicted in Fig. [221 The real-space wave- 
fronts shown in Fig. [2Hlare still Cherenkov-like, since the 
separation between the two branches is relatively small. 

The situation of Fig. [2] and [23 is instead different, 
with a full gap opened between the two LP branches. 
The k-space emission is then concentrated around the 
point k a where the Bogoliubov dispersion touches the line 
lu = u)p. Correspondingly, the near-field pattern shows a 
localized perturbation around the defect, with a peculiar 
stripe pattern. This pattern is due to the interference 
of the Rayleigh scattering at k a and the pump beam at 
k p , so that the wavevector corresponding to the fringes 
is equal to k p — k a . 

VII. CONCLUSIONS 

In conclusion, in this Festschrift paper, we have pre- 
sented a comprehensive analysis of the exotic collec- 
tive excitations of a moving polariton fluid driven by a 



continuous-wave optical pump, which were recently pre- 
dicted in a short Letter [23. We have analyzed in detail 
the interplay between the non-trivial dispersions of the 
polariton Bogoliubov excitations and the onset of single- 
mode (Kerr) or multi-mode (parametric) instabilities. 
We have studied the propagation of the polariton fluid 
in presence of static perturbation potential acting both 
on the photonic and exciton component of the polariton 
field. In particular, we have shown the strict connection 
between the dispersion of the elementary excitations in 
a quantum fluid of microcavity polaritons and the in- 
tensity and shape of the resonant Rayleigh scattering on 
defects. We have pointed out some experimentally acces- 
sible consequences of polaritonic superfluidity for realis- 
tic microcavity parameters. In addition, we have shown 
that in the present non-equilibrium system, it is possible 
to have spectra of collective excitations, which are not ac- 
cessible in systems close to the thermal equilibrium, such 
as superfluid helium or the ultracold atomic condensates. 
These spectra of excitations can be dramatically changed 
by tuning the pump excitation parameters, namely its 
frequency, incidence angle and intensity. We have shown 
a very rich phenomenology in the far-field and near-field 
images of the resonant Rayleigh scattering emission. 
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